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FLIGHT INVESTIGATION OF A SUPERSONIC CANARD 
MISSILE EQUIPPED WITH AN AUXILIARY DAMPING-IN-PITCH 
CONTROL SYSTEM 


By Martin T. Moul 
SUMMARY 


A 609 delta-wing canard missile equipped with а rate-damping system 
to dampen longitudinal oscillations has been flight-tested through а Mach 
number range of 0.70 to 1.85 with a static margin at supersonic speeds 
of approximately 50 percent of the mean aerodynamic chord. Experimental 
transient responses to step deflections of the canard surfaces are pre- 
sented to illustrate the auxiliary damping in pitch provided by the rate 
system. A comparison is made of experimental transient responses and 
theoretical transient responses computed from the linearized equations of 
motion by application of operational-calculus methods and servomechanism 
theory. Some missile aerodynamic characteristics are presented as func- 
tions of Mach number. 


The rate-gyro—servo system, acting through wing-tip controls to pro- 
vide auxiliary damping in pitch, was effective in damping the transient 
responses and caused the oscillations to be almost completely damped at 


the end of 15 cycles. The comparison between measured and computed | 


responses showed that the theoretical method provides an accurate means 
of predicting missile response if experimental stability derivatives and 
rate-system characteristics are used. The experimental data indicated 
that the addition of tip controls to the canard missile caused a varying 
lift-curve slope during the oscillatory response, reduced the pitching 
effectiveness at transonic speeds, but had no appreciable effect on trim 
drag. 


INTRODUCTION 


Investigations of the performance characteristics of an automatically 
controlled canard missile configuration being conducted by the Langley 
Pilotless Aircraft Research Division have indicated that the aerodynamic 
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damping in pitch of the configuration determined from rocket-powered 
model tests (refs. 1 and 2) is inadequate for certain guidance problems. 
A theoretical study (ref. 3) has shown that the missile dynamic perforn- 
ance characteristics are improved by the addition of auxiliary damping 
introduced by а rate-gyro—servo control system through either canard 
fins or wing-tip controls. The present investigation, with the rate- 
demping system acting through wing-tip controls, was undertaken to verify 
experimentally the results indicated in reference 3. 


Auxiliary-damping data are presented as transient responses to a 
step-function control-surface input for a Mach number range of 0.70 
to 1.85. Some longitudinal aerodynamic characteristics are also pre- 
sented for the Mach number range of the flight. 


SYMBOLS 

с wing chord, ft й | = 
с wing mean aerodynamic chord, ft 
Б total wing area in one plane including body intercept, ва ft 
t time, sec 
h altitude, ft i 
" mass, slugs и 5 
W weight, lb р Е 
Іу moment of inertia about Y-axis, slugs-ft? . T 
g acceleration due to gravity, ft/sec? 
V velocity of model, ft/sec | | 
Ус Speed of sound in air, ft/sec 1 | 
М Mach number, V/V. 
q dynamic pressure, 1Ъ/за ft 
Ө angle of pitch, deg | ын "Ж = 
‘ef angle of attack, deg : | 

FREE TAL. 
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Subscripts: 
trim - 

8 

+ 


е 


а, 


с,6,6,а 


control-surface deflection, deg 
pitching velocity, deg/sec 


normal accelerometer reading 


longitudinal-accelerometer reading, deceleration positive 


а 
lift coefficient, Є сов о - зе віп а 5 


a a 
drag coefficient 21 cos а + 20 віп а)-Ї 
і ‚ (2 g 45 


pitching-moment coefficient, 
Pitching moment &bout center of gravity 
а5с 


Laplace transform variable corresponding to the differential 


operator D= а 


dt 


frequency of oscillations, radians/sec 


trim or steady-state condition 
canard control surf&ce 
wing-tip control surface г 


equivalent 


= 


2 


eR 
< 


Чо! 


partial derivatives 
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APPARATUS AND METHODS 


Model Description = 


The missile fuselage consisted of a {-inch-diameter cylindrical sec- 
tion and ogival nose and tail cones. Triangular wings were mounted оп 
the body in a cruciform arrangement. Smaller triangular fins were mounted 
near the nose in line with the horizontal wing as &ll-movable control sur-. 
faces. Тһе model arrangement and principal dimensions are presented in 
figure 1. = 


The triangular wings were swept back 599 32' and had а modified А 
double-wedge airfoil section. The smell triangular Canard fins were swept 
back 60° and had the same airfoil section as the wings.  Half-delta fins 
with the leading edge swept back 59930! and of а double-wedge section were 
mounted at the tips of the horizontal wings as all-movable controls. Тһе 
Wings were constructed of solid magnesium and the two sets of control sur- 
faces of steel. Details of the wing and control surfaces are presented 
in figure 2. It may be noted in figure 2 that the tip controls are slightly 
smalier than the replaced wing tip. A photograph of the model is pre- 
sented as figure 3 and physical characteristics are presented in table I. 


Control Systems 


The auxiliary-damping system consisted of a rate gyroscope, propor- 
tional low-lag pneumatic servomotor, slide valve, air accumulator, reg- 
ulator, purifier, and linkages. All the components, with the exception 
of the accumulator, regulator, and purifier, were located immediately 
rearward of the wings and occupied an -inch section of the fuselage. А 
photograph of the system installation is presented as figure (Қа). 


Rate-gyro signals were transmitted by mechanical linkages as shown 
by the schematic diagram in figure h(b), through the air valve, servo- 
motor, and torque rod to produce deflections of the wing-tip control sur- 
faces proportional to pitching velocity. The motion of the gyro gimbal 
was spring-restrained and damped by two dash pots. A Lanchester damper 
was used to damp the servo motion. The mechanical linkage between the 
servomotor and torque rod was designed to limit the tip controls to max- 
imum deflections of +10°. А block diagram of the missile and auxiliary- 
damping system is presented in figure 4(c). 


The canard fins were deflected in a continuous square wave with 
deflections of 5.2° and -4.8° and a dwell time of 0.75 second, to pro- 
vide pitch disturbances. A hydraulic pulsing ауре вир Ss from ап 
accumulator actuated the controls. 
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Instrumentation 


The model was equipped with an NACA nine-channel telemeter which 
transmitted continuous records of normal (two ranges), longitudinal, and 
transverse accelerations, angle of attack, canard control deflection, 
wing-tip control deflection, total pressure, and a calibrated static 
pressure. A free-floating vane extending from the nose on a sting meas- 
ured angle of attack, and a tube extending below the body measured total 
pressure. The positions of these two instruments аге shown in figure 1. 
The measured angles of attack were corrected to the missile center of 
gravity by the method of reference Ч. 


The model trajectory was determined by a modified SCR 584 type radar 
tracking unit. A radiosonde released at the time of flight measured tem- 
peratures and atmospheric pressures through the altitude range traversed 
by the model. Model velocity was obtained from a CW Doppler velocimeter 
and from total and radiosonde static pressures. 


Launching 


The model was boosted to supersonic speeds by two solid propellant 
rocket motors of 20,000-pounds-seconds total impulse and 3-second burning 
time which were ignited simultaneously. The method of launching was віп- 
ilar to that described in reference 1, with the launching angle being 
changed to 60°. 


Method of Analysis 


In reference 3 several methods were presented for studying missile 
automatic-control-system performance by using frequency-response and 
transient-response characteristics. The method for calculating transient 
responses from the linearized equations of motion is used herein to pro- 
vide a comparison of experimental and theoretical results. The method, 
as adapted for solving by the Reeves Electronic Analog Computer (REAC) is 
presented here. 


The block diagram representing the system is similar to that of 


reference 3, except that the attitude-autopilot loop is excluded and an 
airframe transfer function a/@ is added: 


Rate Control 
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Airframe transfer functions were calculated by using the method of 


operational calculus from the longitudinal equations of motion for two 
degrees of freedom as used in reference. 3. The simplified (9 = © 
в 


Open-loop transfer functions are 


(5838 398 б.) : БЕ | | 


CS басе 
gran 395 ADT. Б dos 345 Uma рү * 57.308 Omg, BY 
шу _ = | 
"lo, 57. xis - Чин Uma + Сота gy 


а = р (2 (ey 
+ Сү о | 
la 51,346 _ 


Experimental stability derivatives of the configuration, reported 
in reference 2, апа flight conditions &nà model characteristics from 
this investigation were used in determining the transfer functions. Val- 
ues of the aerodynamic quantities used in the calculations аге listed 
in table II. Тһе aerodynamic derivatives for a Mach number of 1.81 were 
obtained by extrapolating the data of reference 2. Tip-control data 
were obtained from the results of two tests reported in references 5 
and 


The rate-gyro—servo system was originally а single-degree-of- 
freedom system with the natural frequency changed from 88 to 221 radians 
per second. Because the force provided by the servomotor was only 
30 pounds, the maximum allowable aerodynamic hinge moment was small. As 
в result an inertia damper was added to the servomotor to reduce the mag- 
nitudes of the initial overshoot of the servo response. This modifica- 
tion enabled а higher static sensitivity 54/6 to be used. From а con- 
sideration of hinge moments and pitching velocities, а static sensitiv- 
ity 9t/Ó оғ 0.2 was chosen as the maximum feasible value. This value 
is in line with the results of reference 3 that high static sensitiv- 
ities yield more satisfactory responses than small sensitivities for this 
range of gyro-servo dynamics. m 


With this sensitivity the rate-gyro— servo system was then tested 
and a transient response to а step deflection of the gyro gimbal 


OSEE 


PEN 
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was obtained. The modes of motion were recognized in the response, and 
an equation was fitted to the curve by a trial-and-error method. Тһе 
system transfer function, determined from this equation by the methods 
of operational calculus, is 


St _ _-0.1р(р® - 6,386р - 2,778,600 ) (3) 
Ө (р3 + TOp? + 49,825p + 1,945,000) 


The tip-control deflection 6, was fed back into the system as ап 
5, | 

. The operations in 
E 


solving the equations to determine transient responses to а 6, input 
of 5° were performed by the REAC and the results were obtained from 
recording elements. 


equivalent canard deflection Кӛ, where К = 


Accuracy 
Inaccuracies in the experimental data result from errors in telem- 


eter equipment and radar. The measured quantities are believed to be 
accurate within the following limits: 


Limit of Accuracy 


+0.6 +0.029 +0.01} 
+.6 + , 00} +.002 


The REAC solutions of a and 65+ are believed to be accurate 
within 20.19 and 10.29, respectively. 


RESULTS AND DISCUSSION 


Damping Characteristics 


Effect of auxilia in pitch on model response.- In order 
to show how this rate-gyromservo control system damped the model oscil- 
lations, some experimental angle-of-attack transient responses to a 
5O gtep input of the canard surf&ces &re presented in figure 5 for Mach 


— al 


8 COMES DENTES МАСА RM L52K14b 


numbers 1.81, 1.32, 0.97, 0.71, and 0.72. Résults of REAC calculations 
are also presented to show the model angle-of-attack transient responses 
which could be expected if no auxiliary damping were provided. The 
increase in damping is apparent with the experimental responses being 


almost completely damped at the end of 15 cycles at all Mach numbers. In 


addition to providing а more satisfactory response to autopilot сошшапа 
Signals, the missile with auxiliary damping experienced & reduction in 
maximum normal loads of about 33 percent, an important consideration for 
maneuvering flight at high Mach number. An effect of altitude was obtained 
when the model accelerated to a Mach number of 0.72 during descent. At 

the lower altitude the oscillations were damped in a shorter time as the 
result of a higher model natural frequency. 


Comparison Of experimental and theoretical transient гевропвев.- 


In figure 6 measured angle-of-attack and tip-control transient responses | 
to а 59 step input of the canard surfaces are compared with REAC solutions 
of closed-loop transient responses at Mach numbers of 1.81 and 0.71. "At 
а Mach number of 1.81 the experimental transient response has a higher 
steady-state angle of attack and lower natural frequency than the com- 


puted response; thus a possible difference in Cy, is indicated. Ехрег- 


imental data were not available above a Mach number of 1.45 in refer- 
ence 2, so that errors may have been incurred by extrapolations of the 
derivatives for the calculations to a Mach number of 1.81. The damping 
of the responses is in good agreement indicating the validity of the 
theoretical method for predicting damping characteristics of а missile _ 
equipped with a rate-damping control system. 


The tip-control deflections at M = 1.01 аге in good agreement 
except for the maximum values. The tip controls of the model were limited . 
to deflections of +10° because of hinge-moment considerations. 


At a Mach number of 0.71 the experimental and computed responses are 
in good agreement, except for a small difference in natural frequencies. 


Rate-gyro—servo response characteristics.- The gyro-servo transient 


response presented in figure 7 was obtained by measuring the tip-control 
deflection when a step input was applied to the gyro gimbal. The curve 

was fitted by an analytical expression from which the transfer function 

was obtained. А comparison with the gyro-servo E given 1n figure 11 


А damping mode and increased natural frequency may be noted: 


Because of servo-energy hinge-moment considerations discussed рге-. 
viously, the dynamics of the gyro-servo system of this test were altered 
во as to decrease the maximum control deflections. This difference in 
control deflections may be noted in figure 8, in which model angle-of-attack 
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and tip-control-deflection transient responses to а 50 step input of the 
canard controls аге presented for а Mach number of 1.81. Тһе angle-of- 

attack responses indicate that the auxiliary-damping system of this test 
did not damp the oscillations as well as the original gyro-servo system 

might have with the use of higher initial tip-control deflections. 


Aerodynamic Characteristics 


All the experimental data of this investigation occurred at a reduced 


frequency Шс 5 0.013 апа should be free of unsteady-lift effects. 


Lift-curve slope.- Model lift-curve slope CT. was determined from 


the lift-coefficient апа angle-of-attack time histories by the method of 
least squares and is presented in figure 9. А large difference is noted 
between the data for increasing and decreasing a. This difference may 
be due to a nonlinear lift curve or possibly to a tip-control effect. 
The existence of a nonlinear lift curve was noted in reference 1 and the 
data for the lower angle-of-attack range are presented. This difference 
in HR may have resulted from а tip-control effect since а majority of 


the control responses occur during the times of increasing а. Computa- 
tions were made to determine the magnitude of lift due to rate of con- 
trol deflection but the contribution.of this term to the model lift was 
found to be negligible. Тһе UE at decreasing о is in fair agree- 


ment with that of reference 2. 


Steady-state angles of attack, lift coefficient, and tip-control 
deflection.- The steady-state angles of attack, lift coefficients, and 
tip-control deflections are presented in figures 10, 11, and 12, respect- 
ively, ав a function ОҒ Mach number for the two canard deflections. The 
data show а gradual variation of angle of attack and lift coefficient 
with Mach number with the exception of Mach number: 0.92, where a larger 
trim angle of attack and lift coefficient are noted.  Tip-control deflec- 
tions increased with increasing Mach number from 0.8° at M = 0.72 to 3.50 
at М = 1.86. At the lower speeds there is а considerable scatter of the 
control-deflection data, but the scatter is within the accuracy of the 
data (40.69). 


Steady-state angle of attack per unit control deflection.- From 


the Qirim Curves, the steady-state angle of attack per unit control 
deflection was determined and is presented in figure 13. А comparison 
with erim from reference 2 (for the same center-of-gravity location) 


5 
Абв аль 
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shows that auxiliary damping has little effect upon control effectiveness 
‘at supersonic speeds. “At transonic speeds there was a large reduction 


in “trim with the addition of auxiliary damping. 
8 | 


At four Mach numbers calculated values of “trim from the experi- 
8 


mental data of reference 2 аге modified to include the pitching-moment 
contribution of the tip controls. The test results were larger at super- 
sonic speeds and lower at transonic speeds than the calculated values. 
The indication is that the two configurations differed somewhat in 


either Cm, or Ош, , ав was suggested in the discussion of figure 6. 


Trim drag coefficient.- Trim drag coefficient is presented in fig- 


ure 14 with the trim drag coefficient for the model of reference 2. Тһе 
data show that the addition of rate damping incurred no drag penalty. The 
differences in the values at subsonic speeds may be mainly attributed to 
inaccuracies of Mach number апа the longitudinal accelerometer at low 
speeds. 


CONCLUSIONS 
The results of a flight investigation of & canard missile having a 


static margin at supersonic speeds of approximately 50 percent С and 


equipped with a rate-gyro—servo control ‘system to supplement the inher- 
ent aerodynamic damping in pitch indicate the following conclusions: © 


1. Whereas the canard missile was normally lightly damped in the 
pitching mode, the addition of the rate-damping system to the missile 
caused the transient responses to be nearly completely damped at the end 


of 1} cycles at all Mach numbers. i і З 


2. The experimental and calculated responses were in good agreement 


indicating that the theoretical method provides an accurate means of рге- 


dicting missile response when а rate-gyro—servo control system is used 
to provide additional damping in pitch. 


3. The experimental data indicated a difference in lift-curve slope 
during the oscillatory response which may be attributed to either а non- 
Linear or tip-control effect. Pitching effectiveness as measured by the 
steady-state angle of attack produced per unit control deflection was 


званим L 
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decreased at transonic speeds by the inclusion of the rate-damping control 
system. The pitching effectiveness at supersonic speeds was not affected. 
No differences in trim drag were noted due to the addition of the rate- 
damping control damping. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 


Langley Field, Va. 
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TABLE I 
MISSILE PHYSICAL CHARACTERISTICS 


Weight and balance: 


Weight, lb : 818 авы | 110 

Center of gravity, station in inches ЭГЭЭ ЛЖ” 59.9 

Pitch moment of inertia, slug-ft? .... dx. 27,6 
Wing: 

Area, in one plane including body се ва ft ЭРЭЛ”... 

Mean aerodynamic chord, ft .. Lu de x ri^ m 4... 1,46 

Thickness ratio st winecbody Juncture- "^ "n СЭР” 0.03 
Canard control: 

Exposed area, two fins, sq ft ..... ЭР” . . 0.19 

Root chord, ft : EE ub ws E їз (cw cw. 10:50 

Thickness ratio at fin-body juncture "A .. M 0.03 
Tip control: 

Exposed area, two d ЕСЕП ж зок ко вою „ ve. ЧЭ... 

Root chord, ft . . * е т е » е е е ч а е ь е “ е е . . е . 0,53 

Thickness ratio Ж. бй ыы ы. AO EL S жом UE. Er ue es cee s DER ЭГ 0.03 

еее Siati кб 
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TABLE II 


AERODYNAMIC DERIVATIVES FOR THEORETICAL RESPONSES 


[The aerodynamic derivatives, as used in the calculations 
and tabulated below, are per degree measure. The static 
margin was approximately 50 percent 6 at supersonic 


speeds. 


Mach С 
шы ы|=®|® | е | | | №6 


т 


Тір control шин?! 


Canard control 59° 52 i 
м миний D 
Pal 


? Ер) zi s ЧЕР БЭ Б T 
— n гү ч = Ш ! iR 
^ ! 
Plan view aj a 
- Station 
о 52 „В 59.9 88.5 114.2 
Angie-of-attack Indicator $ 
31.0 
= — --- 5 —— ®———— === - Вв" 


Side view 


Total~pressure pickup 


Figure 1.- Model arrangement. А11 dimensions аге in inches. 
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Wing 


Figure 2.- Detail of wing апа control вигРасев. All dimensions are in 
inches. 
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Figure 3.- Missile plan view. 
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(а) Rate-gyro-servo system installed in missile. 


Figure 4.- Rate-gyro-servo system. 
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Air supply 
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Air valve ау! 
h 


JK mas servomotor 37 


Inertia damper i 
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(b) Schematic diagram of components and linkages. 


Rate control 


(с) Block diagram of auxiliary-damping system incorporated in missile. 


Figure 4.- Concluded. 
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а, deg 


«Р 


- 4 


с, deg 


Figure 5.- Effect of auxiliary damping in pitch on model response. 
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(a) М-1.81; h = 5,000 feet. 


ЩЩ 
ЕЕРЕЕ 
ЕА ТТТ 
АТТ 
HAN eS 


/ 
А ГЕК 
тү ТҮ 
тг ЇГ 
НН 
ДЕЕЕЕЕЕЕЕШ-ЕЕНн; 
А95 


Time,t,sec 


(b) М = 1.32; h = 9,500 feet. 
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(с) М = 0.97; h = 13,000 feet. 
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Time,t,sec 
(а) М = 0.71; В = 17,000 feet. ` Й 


Figure 5.- Continued. 
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Test (auxiliary damping) 
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(e) М = 0.72; h = 6,000 feet. 


Figure 5.- Concluded. 
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ec, deg 


Time,t,sec - = 


(а) M = 1.81. B | Е 


Figure 6.- Measured апа computed transient responses with rate damping. 
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Figure 6.- Concluded. 
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Figure 7.- Experimental rate-gyro-servo transient response to а unit 


pitching-velocity input. 
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Figure 8.- Effect of rate-system dynamics on calculated model responses 
at a Mach number of 1.01. 
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Figure 10.- Variation of trim angle of attack with Mach number. 
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Figure 11.- Variation о? trim lift coefficient with Mach number. 
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Figure 12.- Variation of trim tip-control deflection with Mach number. 
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Figure 13.- Variation of trim angle of attack per unit control deflection 
with Mach number. 
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Figure 14.- Variation of trim drag coefficient with Mach number; 
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